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A geometric potential Vc depending on the mean and Gaussian curvatures of a surface E arises 
when confining a particle initially in a three-dimensional space Q onto E when the particle Hamil- 
tonian Hn is taken proportional to the Laplacian L on fi. In this work rather than assume Ha oc L, 
momenta P v Hermitian over Q are constructed and used to derive an alternate Hamiltonian H v . 
The procedure leading to Vc, when performed with H v , is shown to yield Vc = 0. To obtain a 
measure of the difference between the two approaches, numerical results are presented for a toroidal 
, model. 
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1. INTRODUCTION 

The study of quantum mechanics on curved surfaces has been a subject of theoretical effort for decades 0, 0, EJ 0, 
IE. IM 1^ [10, [lH E3, EH EH HM. H^, . The approach as pioneered in [J, H H first defines a surface 
v) and employs a coordinate q to describe excursions at any point normal to X. The Hamiltonian appropriate 
to the three dimensional coordinate system f2 characterized by a metric gij is then defined by Hq = — iV 2 with 
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As will be briefly reviewed below, Eq.(l), along with imposing conservation of the norm in the limiting case where 
the particle is constrained to motion on X, leads to a geometric potential Vc- This paper with is concerned with how 
Vc is modified when instead of defining Hq through Eq.(l), another choice of Hamiltonian, one for which surface and 
1 1 ^ \ normal momenta are required to be Hermitian over Q is adopted. 

The remainder of this paper is organized as follows: in section 2 the geometric potential is derived from Eq.(l) 
for a cylindrically symmetric geometry fl c . This material has appeared other places, but is included here to keep 
this work as self contained as possible and to establish some notation. A symmetric geometry was chosen for ease 
of presentation; the salient points for more general cases remain intact. In section 3 the classical Hamiltonian for 
a particle in Q c is written, Hermitian momenta on Q c are found, and standard quantum mechanical rules employed 
to write the first quantized Hamiltonian. Finally, using the limiting procedure detailed in section 2, the geometric 
potential is re-derived. Section 4 presents numerical calculations for a toroidal system, and section 5 is reserved for 
■ conclusions. 

2. THE GEOMETRIC POTENTIAL; A BRIEF REVIEW. 

The derivation of Vc will be performed for a cylindrically symmetric geometry J7 C with Hq oc L. The extension to 
the general case is straightforward and the central results below will be not be affected by the restriction to azimuthal 
symmetry. 

Let e p ,e^ be cylindrical coordinate system unit vectors defined through 

= cos(/) i + sin0 j (2) 

= — sin0 i + coscj) j (3) 
and let a cylindrically symmetric surface S c be described by 



r(p» = pe p + S(p)k. 



(4) 
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S(p) gives the shape of the surface. Any point near S c may be reached by travelling normal to S c up (or down) a 
distance q via 

x(p, <j>, q) = r(p, <j>) + qe n (5) 

with e„ everywhere normal to the surface. The metric near the surface derived from Eq. (5) is 0] (subscripts here 
indicate differentiation) 



dx 2 = Z* 



qS pp 
z 3 



dp 2 + P 2 
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Z 2 [l + qhYdp 2 + p 2 [l + qk 2 ] 2 d(j) 2 + dq 



(7) 



with Z = yl + S 2 . The Laplacian can be found straightforwardly from Eq. (7) but is not essential to the main 
arguments that follow. 

Now consider a situation where a large confining potential V n (q) everywhere normal to E c acts to restrict the 
particle to the surface. The derivation of the geometric potential is independent of its specific form [2l| , but however 
chosen it causes the particle to sit on £ c . In this limit the wave function decouples into surface and normal parts 



*(P, </>,<?) -* Xa(P, <l>)Xn(<l)- 

Conservation of the norm in the decoupled limit is imposed by insisting upon the condition 

\*\ 2 FdEdq = \ Xs \ 2 \Xn\ 2 cEdq 
with F = 1 + 2qh + q 2 k and dY, the surface measure, h and k are the mean and Gaussian curvatures given by 

h = ~(fci + k 2 ), 



k = kxk%. 



Inserting the expression given by Eq.(8) 



* = 



XsXn 



(8) 
(9) 

(10) 
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into the time independent Schrodinger equation — iV 2 \t = E^>, and taking the q — + limit after performing all q 
differentiations, yields the pair of equations (h = m = 1) 
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2[Z 2 dp 2 Z 2 pdp Z 3 dp <•• 
ld 2 Xn 



Xs = E s Xs 



The curvature potential Vc is 



2 dq 2 



Vc 



+ V n {q)Xn = E n Xn- 



-[h 2 -k]; 



(13) 
(14) 

(15) 



seeing specifically how Vc is generated will prove important in the following section. It arises from the g->0 limit of 
the 



^+^y s (p^) X n(q)F-^ 2 



(16) 



term in the full Laplacian . The Z dependence appearing in the kinetic energy operator given by Eq. (13) occurs 
independently of the q degree of freedom. The effect on the spectra of some simple systems due to metric factors in 
the operator was shown to be small in 0,0] for some parameterizations of S(p), but it is also possible for the effects 
to be large ^ . 
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3. THE HAMILTONIAN WITH HERMITIAN MOMENTA 

In this section a Hamiltonian H n quadratic in momenta P v Hermitian over fi c is developed. The procedure followed 
comprises standard quantum mechanical rules, to wit: 

a. the classical Lagrangian is written 

b. the canonical momenta are determined 

c. the classical momenta are replaced with their first quantized form P v Hermitian on Q c . 
The Lagrangian for a particle in a geometry described by Eq. (6) is 

L = ^(Z 2 [l + k lQ ] 2 p 2 + p 2 [l + k 2 q] 2 ^ 2 + f). (17) 

With P n = dL/dq v , the classical Hamiltonian for this geometry is 

1 / P 2 P 2 \ 

H= 2\Z 2 [llk iq ] 2 + p 2 [l + k 2 q] 2 +P 'J- (18) 

To proceed with first quantization via P n — > P v , it is necessary to produce momentum operators Hermitian on a 
geometry described by gij . The required relation is 

Anticipating the q — > limit, it is useful to note the differentiations done by the surface momenta will pass over any 
appearances of the q variable, making it possible to set q = in both of them without affecting the final result. The 
momenta in this limit are 

M(| + >+<4 (22) 

H v is straightforward to obtain from the three momenta listed above, but it is Eqs. (16) and (22) that will be focused 
upon to generate the central result of this paper. Write 

A2 d 2 2„ ,.9 1 fdh \ 1 dF„ 1 ,, 2 

- P « = W 2 + F {H + qk) 0q + F U + V - ^ + ^ + + qk) ] ^ 
upon taking the q — > limit the result is 

- p ' = w 2+2h -k +k - h2 - (24) 

The key point is simple; Eq.(24) includes the operator that appears in the Laplacian formulation of the problem as 
shown in Eq.(16). It eventually yields Jjp- plus h 2 — k, from which it is trivially seen that the geometric potential 
cancels in the Hermitian formulation. 



4. A NUMERICAL EXAMPLE 

In this section the spectra and wave functions of the Schrodinger equation for a particle on a toroidal surface T 2 
of major radius R and minor radius a will be compared using the two formalisms detailed above. There are good 
reasons to choose such a system, not the least being calculating system observables is tractable. Additionally, toroidal 



4 



nanostructures of varied types have been fabricated so there is some relevance to real devices, a point which will be 
elaborated upon in the conclusions section below. 

The surface of a doughnut shaped torus has may be parameterized by pflj 

x(0, (/), q) = (R + a cos#)e p + a sinflk + qe n . (25) 

Applying d to Eq.(25) gives 

dx = (a + q)d9eg + [R + (a + q) cos9]d(j)e < f > + dqe n (26) 

with eg = —sm9e p + cos#k and e„ = x eg. 

Let a = a I R, (3 = 2Ea 2 ; applying the formalism described in section 2 and making the standard separation ansatz 
for the azimuthal part of the wave function ^(0, <p) — ip(9)exp [ivcfi] yields the Schrodinger equation with Vc present 



d 2 ip a sin dip {v 2 a 2 - \) 



d9 2 [l + a cos 9} 88 [1 + a cos 9} 2 
The method of section 3 gives the Schrodinger equation 



iP + f3iP = 0. (27) 



d 2 ^ a sin 8 dip [ V 2 a 2 + \(a 2 -\)] 1 

-W - [1 + Q cos0] ^ [I + a cos6] 2 ^ + ^-4^ = °' (28) 

Even before consideration of the solutions of Eq. (27) and (28), there is an immediate distinction between the two. 
The cancellation of the azimuthal kinetic energy term at discrete values of a occurs in the former at values of 

a = -L (29) 
2v 

and in the latter at 



The ramifications of Eq.(29) and by implication Eq.(30) have been discussed elsewhere. Here it need be noted that 
while there appears a term in the Hermitian formalism that is identical to Vc derived from the Laplacian method, it 
is an artifact of the high degree of symmetry of the torus and may take a very different form in the general case. 

Numerical solutions of Eqs. (27) and (28) were determined by a basis set expansion in Gram-Schmidt functions on 
T 2 (2j|. Eigenvalues and wave functions are given in tables 1-3 for several low-lying states with a = 1/3, 1/2 and 2/3. 
The results indicate clearly that the two Hamiltonians can differ, particularly at the "magic" a given by Eq. (29) or 
Eq. (30), and when the torus fattens to values of a ~ 2/3. 



5. CONCLUSIONS 



This work demonstrated that a Hamiltonian quadratic in momenta Hermitian on a curved surface and consistent 
with standard quantum mechanical rules leads to a vanishing geometric potential Vc- It is interesting that Golovnev 
has recently employed an alternate (Dirac quantization) procedure for which Vc = on a sphere 26]; however, 
it is not clear to the author how the DQ method could be easily implemented on T 2 . Numerical results were 
presented to show the differences in spectra and wave functions between the two Hamiltonians discussed here can be 
substantial. It is possible to speculate as to whether advances in nanostructure fabrication may perhaps allow for 
experimental determination of which (if either) Hamiltonian studied here is the more correct description of curved 
surface nanophysics. Given that toroidal structures have been fabricated and are "clean" with respect to curvature 
effects (in that there is enough symmetry to make the problem interesting but not so much to make it trivial), it is 
conceivable that existing devices may be employed to settle the very fundamental issue of what is the appropriate 
Hamiltonian for a particle constrained to a surface. 



6. ACKNOWLEDGEMENTS 



The author would like to thank Mr. Lonnie Mott for useful discussions. 



■5 



H. Jensen and H. Koppe, Ann. of Phys. 63, 586 (1971). 
R. C. T. da Costa, Phys. Rev. A 23, 1982 (1981). 

R. C. T. da Costa, Phys. Rev. A 25, 2893 (1982). 

M. Burgess and B. Jensen, Phys. Rev. A 48, 1861 (1993). 

P. Exner and P. Seba, J. Math. Phys. 30, 2574 (1989). 

S. Matusani, J. Phys. Soc. Jap. 61, 55 (1991). 

P. Duclos and P. Exner, Rev. Math. Phys. 7, 73 (1995). 

J. Londergan, J. Carini, and D. Murdock, Binding and scattering in two dimensional systems; applications to quantum 

wires, waveguides, and photonic crystals ( Springer- Verlag, Berlin, 1999). 

J. Goldstone and R. L. Jaffe, Phys. Rev. B 45, 14100 (1991). 

P. Ouyang, V. Mohta, and R. L. Jaffe, Ann. of Phys. 275, 297 (1998). 

I. Popov, Phys. Lett. A 269, 148 (2000). 

S. Midgley and J. Wang, Aus. J. Phys. 53, 77 (2000). 
I. J. Clark and A. J. Bracken, J. Phys. A 29, 4527 (1996). 
P. C. Schuster and R. L. Jaffe, Ann. Phys. 307, 132 (2003). 
M. Encinosa and B. Etemadi, Phys. Rev. A 58, 77 (1998). 
M. Encinosa and B. Etemadi, Physica B 266, 361 (1998). 
K. Lin and R. L. Jaffe, Phys. Rev. B 54, 5757 (1996). 

A. Chaplik and R. H. Blick, New J. Phys. 6, 33 (2004). 
S. Qu and M. Geller, Phys. Rev. B 70, 085414 (2004). 

B. Nilsson, cond-mat/0103029. 

M. Encinosa, L.Mott, and B. Etemadi, Phys. Scr. 72, 13 (2005). 
M. Encinosa, unpublished. 

M. Encinosa and B. Etemadi, Found. Phys. Lett. 16, 403 (2003). 
M. Encinosa and L.Mott, Phys. Rev. A 68, 014102 (2003). 
M. Encinosa, physics/0501161, submitted to App. Math. Comp. 
A. Golovnev, quant -ph/0508044. 



6 



TABLE 1: Eigenvalues and wave functions corresponding to the solutions of Eqs. (27) and (28) for a — | (the (27r)~i 
normalization from from the <f> dependence is omitted.) Terms not shown are at least an order of magnitude smaller than those 
listed. 



p 


a = ± Vc^O 


-.2834 
-.1528 
+.1968 


*o = -4082 - 0.0776cos<9 

*i = [.4049 - 0.0421cos6']e ± ^ 

* 2 = [-4060 - 0.0525cos6»]c ±2i< ^ 


P 


a = i Vc = 




+.1267 
+.4735 


*o = -.4075 - 0.0686cos<9 
*i = [.4038 - O.OSSScos^e* 1 * 
* 2 = [-3869 - 0.0593cos(9]e ±2i< ^ 



TABLE 2: As per table 1 with a=\. The middle state is qualitatively very different due to a = \ defining a "magic radius" 
when Vc = 0. 



p 


a = i V c ± 


-.3511 


*o = -4230 - 0.1470cos<9 





*i = e ±10 


+.6288 


* 2 = [-.3293 - O.1935cos(9]e ±2i0 


P 


a = ± Vc = 





*o = -4191 - 0.1072cos<9 


+.3192 


*i = [.3912 + O.O293cos0]e ±i0 


+.9212 


* 2 = [-.3289 - O^OOScos^e* 210 



TABLE 3: As per table 1 with a = |. The second excited state for the case when Vc = becomes a negative parity function 
in 9. 



p 


a = i V c 


-.5947 

+.2024 

+.4377 


* = 0.4597 - 0.3406cos<9 - O.1000cos2<9 
*i = [0.3570 + O.lO79cos6>]e ±i0 
* 2 = 0.5930cos(9 - 0.1363cos2(9 


P 


a = I Vc = 


3.4(HT 4 ) 

+.5397 

+1.0000 


*o = -0.4369 + O.1516cos0 
*i = [0.3405 +O.1446cos6>]e ±i0 
*2 = -0.5888sin(9 - 0.0992sin26» 



